NASA Technical Memorandum 84356 Airborne Astronomy Program 

Medium Altitude Missions Branch 
Preprint Series No. 9 


NASA-TM-84356 

tf 8300 UK? 5 


The Determination 
of Cloud Masses and Dust 
Characteristics from 
Submillimetre Thermal 
Emission 

Roger H. Hildebrand 

LIBRARY COPY 

f/i 4. Y 9 - 1983 


April 1983 



NASA 

National Aeronautics and 
Space Administration 



LANGLEY RESEARCH CENTER 
LIBRARY, NASA 
HAMPTON, VIRGINIA 


NASA Techncia! Memorandum 84356 


Airborne Astronomy Program 
Medium Altitude Missions Branch 
Preprint Series No. 9 


The Determination of Cloud Masses 
and Dust Characteristics from 
Submillimetre Thermal Emission 

Roger H. Hildebrand, Department of Astronomy and Astrophysics, The Enrico Fermi 
Institute, The University of Chicago, Chicago, Illinois 


NASA 

National Aeronautics and 

Space Administration 


Ames Research Center 

Moffett Field. California 94035 


THE DETERMINATION OF CLOUD MASSES AND DUST 
CHARACTERISTICS FROM SUBMILLIMETRE 
THERMAL EMISSION 


Roger H. Hildebrand 

Department of Astronomy and Astrophysics, 
Department of Physics, The Enrico Fermi Institute, 
The University of Chicago 


Subject Headings : interstellar: matter-- i nterste 1 1 ar : 


abundances 


Accepted for publication in Vol. 2k, No. 3 of 
THE QUARTERLY JOURNAL OF THE ROYAL ASTRONOMICAL SOCIETY 




* 


Received: 1982 December 1*4 


The Determination of Cloud Masses and Dust Characteristics from 
Submillimetre Thermal Emission 

Roge.fi H. Hlld.tbfw.nd 

Department of Astronomy and Astropnysi cs , Department of Physics, 

The Enrico Fermi Institute, The University of Chicago 
5630 South Ellis Avenue 
Chicago, Illinois 60637 
U.S.A. 

ABSTRACT 

The purpose of this paper is to review the principles by which the 
dust masses and total masses of interstellar clouds and certain charac- 
teristics of interstellar dust grains can be derived from observations 
of far-infrared and submillimetre thermal emission. To the extent 
possible, the discussion will be independent of particular grain models. 


I INTRODUCTION 

The construction of large infrared telescopes, such as the NASA 
Kuiper Airborne Observatory, the NASA Infrared Telescope Facility, and 
the United Kingdom Infrared Telescope, and the development of sensitive 
far-infrared and submillimetre photometers have made it possible to 
measure thermal spectra even for very cool clouds (T <20 K) . 
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In nearly all cases, the optical depths, t, determined from the 
spectra at wavelengths longer than ^200 ym are found to be very much 
less than one. Hence the observed flux densities sample with equal 
efficiency the emission from all depths in a cloud. It is especially 
that characteristic of submillimetre flux measurements which one seeks 
to exploit in determining the characteristics of clouds and grains. 

The fundamental difficulty in determining the mass of any cool 
cloud which is opaque to starlight is that the principal ingredient, 
molecular hydrogen, is invisible; the kinetic energies of the molecules 
are far below the minimum energy required to produce observable hydrogen 
transitions. The methods commonly used to overcome that difficulty are: 

(a) Kinematics . When gradients in Doppler shifts can reasonably be 
attributed to rotation of gravitationally bound clouds one may 
estimate the mass necessary to explain the observed radial 
velocities. Often, however, the observed gradients can be 
produced by ejection of unbound material from compact regions 
within the clouds. 

(b) Extinction . When the opacity of a cloud is low enough to 
permit observation of background stars, one may estimate the 
column density by assuming that the relationship between ex- 
tinction and hydrogen column density measured in the intercloud 
medium can be applied with suitable precautions to the cloud. 
This method can be extended to greater column densities by 
measuring extinction in the near-infrared (e.g. K band; Z.2 pm) 
and applying measured relationships between near-infrared and 
visible extinction. 
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(c) Molecular line intensities. When the column density, N(x), of 


some molecule, x, is obtained by observing a particular emission 

line, one may estimate the hydrogen column density, N ( H 2 ) > by 

assuming a ratio N(H,>)/N(x) equal to that measured in a thinner 

cloud for which NtH^) has been found by the extinction method. 

This method is valid cily when the excitation temperatures are 

known ard when the optical depth for the observed transition is 

<1. For cold clouds, these conditions are often not satisfied 

for easily observable molecules such as C °0. For example, 

large variations from cloud to cloud in the apparent ratio of 

N ( 1 3c 1 ^0 ) to N(^C^O) can be produced by differences in line 

widths and corresponding differences in optical depths for the 
13 16 

observed CO transitions. The relationship between carbon 
monoxide abundance and visual extinction has been studied in 
detail by Frerking, Langer, and Wilson (1982). 

Cloud masses may also be obtained from estimates of number 
densities using intensities of col 1 isional ly excited lines such 
as the 2 mm and 2 cm lines of ^CO, with kinetic temperatures 
derived from CO line intensities, and velocity gradients deter- 
mined from emission line widths. The cloud volume is estimated 
from the size projected on the sky. This technique requires 
extensive cloud modelling. 

(d) Thermal emission . Thermal emission can be used as a measure 

of dust volume , or, if a grain density is assumed, as a measure 
of dust mass. The computation of the volume from the flux den- 
sity depends on the determination of grain emissivity. Even 
without knowledge of the grain emissivity, one may determine aloud masses 
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from thermal emission if the extinction method is used to 
estimate hydrogen column densities in calibration objects. 

The thermal emission method also has limitations, as will 
be discussed, but it is the only known method for determining 
the dust mass of a visually opaque cloud without assuming a 
ratio for the masses of dust and gas. As a method for deter- 
mining total mass it is independent of the others except 
insofar as extinction may be used to calibrate both this and 
the molecular line methods. 

The difficulty in determining the composition of interstellar dust 
is the absence of positive identifying features in either the extinction 
curve or the thermal emission spectrum. The absorption peaks are broad 
and are influenced not only by chemical composition but also by crystalline 
structure, grain mantles, impuri ti es, and grain size. Great effort has been 
directed toward the discovery of plausible grain models which can account 
for the extinction curve and for the observed scattering and polarization 
of starlight. In judging whether a grain model is plausible, one favors 
constituent elements which are depleted in the gas phase below assumed 
cosmic abundances. One would prefer to avoid such judgements and instead 
to improve the knowledge of cosmic abundances by direct analysis of inter- 
stellar dust, the most abundant form of observable solid material. 

Although no sure analysis has been devised, the range of acceptable 
grain models can be reduced by measurements of thermal radiation. Emission 
spectra can provide mean values of emissivities , and when combined with 
extinction curves, mean values of dielectric constants. Moreover, the 
polarization spectrum of emitted radiation may provide tests for the dielec- 
tric and magnetic properties of individual grain species. 

We begin with a discussion of dust and cloud masses (sections 2 and 3) 
followed by a discussion of grain properties (dielectric constants, 

5 


emissivities, and composition; sections 4, 5 and 6), and finally by a 
table of numerical values of submillimetre emission parameters (section 7). 


2 DUST MASS 

To the precision now available, one finds emperically that the sub- 
millimetre spectrum of a cool cloud can be fit by the product of the black- 

o 

body intensity, B(v,T) and an emissivity Q(v) which varies as v , where 
g increases from ^1 at v % 50 cm -1 (A 200 pm) to £2 by v % 10 cm" 1 
('vlOOO pm), (e.g. Erickson, et al_. 1981; Schwartz 1982). The spectrum 
thus gives a measure of a nominal temperature, T, falling within the range 
of temperatures of the individual dust grains. The determination of the 
dust mass of a cloud at a known distance is based on measurements of the 
flux density and T, and of Q(v) and the size of the grains. 

We consider first an idealized cloud composed of spherical grains of 
uniform size, composition, and temperature. The cloud is optically thick 
to starlight, is optically thin to submillimetre radiation, and is heated 
by internal processes. We shall later consider the effects of changes in 
these grain and cloud characteristics. 

2.1 Idealized ctoud . The flux density, F(v), from a cloud at a distance 
D containing N spherical dust grains each of cross section o, temperature 
T, and emissivity Q(v), is given by 
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( 1 ) 


F(v) = N[a/D 2 ] Q(v) B(v,T). 


The volume of dust in the cloud is given by 


V = Nv 


( 2 ) 


where v is the volume of an individual grain. Eliminating N from these 
equations, one obtains 


If one assumes a grain density, p, one obtains an expression for the 
dust mass = V p , or 


where a is the grain radius (Hi 1 debrand et aK 1977). 

It may appear that this expression is applicable only for clouds 

with uniform spherical grains and that the radius, a, must be known for the 
grains in each cloud under consideration. We show next that such is not 
the case. (We continue, until section 2.5, to assume uniform composition 
and temperature.) 

2.2 Gruzin ilzo. . The extinction curve for starlight passing through 
interstellar material is not consistent with attentuation by grains of 
uniform size (e.g. Mathis and Wallenhorst 1981). Clearly the quantities 
a and Q = Q(a,v) inequation (4) must be suitably chosen averages. 

The problem would appear to be complicated by the increase in a, and 
hence a change in Q(a,v), in dense clouds as indicated by peculiar ex- 
tinction measurements (Carrasco, Strom and Strom 1973; Coyne, Gehrels, 


V = [F(u)D 2 /B(v,T)][v/a]/Q(v). 


(3) 


M d = [F(v)D 2 /B(v,T)][(4/3)«/Q(v)]p 


(4) 
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and Serkowski 1974; Serkowski , Mathewson and Ford 1975). Fortunately, 
however, for \ >> a, a simplification is imposed by the Kramers-Kroni g rela- 
tions which can be used to connect Q with a and with the complex dielec- 
tric constant, e(v), of the grains at a given frequency v (Jackson 1961; 
Andriesse 1977). For spherical grains, one may write the relationship 
in the form 


Q(v) = 8ira(v/c) Im {-[e(v) - l]/[e(v) +2]}. (5) 

The quantity a/ Q(v) is thus independent of a for a given value of e(v). 
Hence, if values a and Q are correctly determined in some medium, the 
corresponding ratio a/ Q(v) can be used with equation (4) to obtain the 
dust mass of any cloud with the same value of e(v). Moreover, if a/Q(v) 
is independent of a, the flux density from a grain must be proportional to 
oa ( cf . equation 1); that is, to the grain volume., v. 

Although these conclusions drawn from equation 5 are strictly true 
only for spherical grains of uniform composition, they will be important in 
our consideration of non-spherical grains and grains in which ice mantles 
surround refractory cores. 

When an "average" 0(v) is measured as described in section 5, the 

value for each grain is weighted by the contribution of the grain to the 

flux density and hence by the grain volume, v (or a ). The corresponding 

value of a is weiqhted in the same manner; for a size distribution n(a), 

3 3 

the weighted average is fan(a)a da/fr\(a)a da where the integrals extend 
over all values of a. 

Several distributions have been proposed. Mathis, Rumpl , and 
Nordsieck (1977) obtain a fit to the extinction curve in the visible portion 
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in the range 


of the spectrum using a distribution n(a) « a 
0.01 ym _< a _< 0.25 gm. Honq and Greenberg (1980) propose a core-mantle 
grain model with a distribution favoring somewhat larqer radii. From 
ultraviolet scattering data, Witt (1979) infers a distribution including 
smaller radii. Using the Mathis, Rumpl , Nordsieck model, we obtain a 
weighted average radius a = 0.1 ym. 

Ideally one should obtain a and Q(v) in the same cloud. For 
the present we must be content with an a from the intercloud medium and 
a Q(v) from a reflection nebula (section 5). 

From the arguments we have advanced, we conclude that in otherwise 
ideal clouds, variations in a from grain to grain and from cloud to 
cloud do not influence dust mass determinations if corresponding values of 
a and Q(v) have once been established. 


2.3 Grain sha£e_ . The polarization of starlight indicates the presence 
of aligned non-spherical grains. One must therefore consider the extent 
to which the observed flux density from a fixed volume of dust depends on 
the shape of the grains and the extent to which the substitution of (4/3 )a 
(equation 4) for v/a (equation 3) depends on shape where a is deter- 
mined from the extinction curve assuming spherical grains. 

The first question is already answered by equation 5 and the conse- 
quent dependence of F(v) on grain volume if one adds that the equation 
is very nearly true even for greatly elongated or flattened shapes (Purcell, 
1969). Unless one suspects grains in which the ratio of the axes is greater 
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than about five one may safely ignore the effect of shape on flux density 
until Q can be measured with ten times better precision than has yet 
been achieved. 

The second question, the use of an a from the extinction curve to 
find v/o for non-spherical grains, has been answered by Greenberg and Hong 
( 1975) who note "that if one were to take a set of spherical particles and 
elongate them, the shape of the wavelength dependence of extinction by the 
randomly oriented elongated particles would be qualitatively like that of 
smaller spheres. Thus, in order to obtain the same wavelength dependence 
of extinction, we should have started with larger spheres (but smaller numbers) 
in order to maintain the same total extinction. Thus, the total volume of 
material [derived for transparent clouds] turns out to be essentially shape- 
invariant because it is constrained not only by the amount of visual extinc- 
tion but also by the wavelength dependence." 

Aside from .consideration of techniques for estimating grain size if 
one assumes or somehow measures a mean projected cross section, o, of a 
grain, one may still ask how the ratio v/o in equation 3 depends on shape. 

The answer to this purely geometrical question is readily found for any as- 
sumed convex shape using the relationship between mean projected area and 
total surface area derived by Schwartzschi Id (see Appendix). The ratio v/o 
remains within a factor two of its maximum value, (4/3)a, for all regular 
polyhedrons and all circular cylinders of length between ^1/5 and ^7 times 
the diameter. (Here, we define a to be [o/n] ? ). 

We conclude that unless the bulk of dust is made up of very thin 
needles or flakes, the shape of the grains will have little influence 
on the determination of dust volume. 


10 


2.4 Gruzin mantle ^ . In tenuous clouds exposed to starlight one expects 


grains to be composed entirely of refractory materials. In cool dense 
clouds, an absorption feature at 3.1 pm gives evidence for the formation of 
ice mantles where the mass of ice is usually found to be small compared to 
the mass of the refractory cores (Merrill, Russell, and Soifer 1976). 

The effects of chanoes in surface material have been examined in detail 
by Aannestad (1975). The effect on extinction in the near and mid-infrared 
can be very great, but the effect on thermal emission at submillimetre 
wavelengths is relatively small: as we have seen (section 2.2), the entire 

volume of a grain contributes to the emitted flux. The dust mass derived 
from equation (4) should not be significantly changed by mantles unless 
the mantles make up a major portion of the volume. 


2.5 GsuUn mtxXuAej > . The temperature of a grain may depend significantly 
on its composition. Clouds made up of several grain species will exhibit 
composite thermal spectra with nominal temperatures, T, dependina on 
abundances and temperatures for the individual species. For grains of 
widely differing emissivi ties , the temperatures may differ by almost a 
factor of two (Mezger, Mathis, and Panagia 1982). 

If the grains of different species also differ in size, then the 
value of a for the mixture may not be appropriate for grains emitting, 
say, on the low frequency side of the spectrum. This does not mean that 
one must know a frequency dependent value of a: the empirical relationship 

Q(v)*v e allows for the frequency dependence of any parameter affecting Q. 
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What it does mean is that the "average" value of Q should properly be 
determined at a frequency, v , such that the various species contribute 
to the emitted flux in proportion to their volumes; that is, such that the 
product B(v,T)Q(v)/a is the same for each species. Such a frequency will 
exist for two-component mixtures. The best one can do without assuming 
a model for grain composition is to determine Q somewhere near the center 
of the emission spectrum. 

A further complication for mixtures of grains has to do with exposure 
to ultraviolet radiation. Since the ratio of ultraviolet absorptivitv to 
infrared emissivity may differ from one species to the next, the distribution 
of temperatures among the various species may be different in clouds having 
the same mixture of grains and the same mean temperature but different 
exposure to ultraviolet sources. Hence the mean value of a/Q determined from 
submillimetre emission by grains in an optical reflection nebula may riot 
be accurate for a dark cloud of the same composition where the value (a/Q).. 
for the ith species is weighted differently by the corresponding values of 
B(v,T.j). Ultraviolet exposure may also increase the weighting of very 
small grains (a ^ 0.01 mm), which may undergo large temperature fluctuations 
when absorbing single energetic photons (Purcell 1976; Aannestad and 

Kenyon 1979). 

The factors which alter grain temperatures including fluctuations 
of individual grains, differences in the emissivities of the various grain 
species, differences in ultraviolet exposure, and large scale gradients 
due to differences in exposure to heat sources, all tend to broaden the 
observed thermal peaks. With measurements of the shapes of the peaks for 
selected objects one might place upper limits on some of these factors. 
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Even with such limits* the problems related to mixtures of grains 
with different emissivities are likely to remain more troublesome in prin- 
ciple than those related to variations only in size or shape or to the 
growth of mantles; but since the value of Q is still known at best to a fac- 
tor of three for any cloud in any environment, it is likely that mass deter- 
minations can be improved by better measurements of Q whether or not the 
measured Q nicely corresponds to the average of the species according to 
their volumes. 


2.6 ExteAymLly [itaXzd cloucfa . In the special case of clouds heated only 
from the outside by short wavelength radiation for which x >> 1 , the far 
infrared brightness is determined by the incident radiation and the surface 
properties of the cloud; not by the column density of the dust. For example, 
Keene (1981) finds approximately the same surface brightness for many 
globules despite variations in size. We shall leave the analysis of 
externally heated clouds to papers dealing especially with that problem 
and with energy transport in clouds (e.g. Leung 1975; Mezger, Mathis and 
Panagia 1982). 
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2.7 Conclusion to discussion o£ dust mss_ . In clouds which are opaque 


to starlight, the only known method for determining dust mass without 
assuming a ratio of gas to dust is to measure thermal emission. The 
method depends on ar assumed grain density and on an estimate of the distri- 
bution in grain size (in a reference object) based on the extinction curve, 
but is otherwise insensitive to grain models. Changes in the size and shape 
of the grains and the growth of mantles are unimportant except in extreme cases. 
The accuracy of the dust mass determinations is now limited primarilyby the un- 
certainties in grain size and emi ssivi ty. In particular cases, the accuracy may 
be 1 imi ted by large temperature variations along the line of sight: the possi- 

bility of underestimating the masses of cool components tends to increase 
with increases in the nominal temperature derived from the peak of the 
spectrum. 

As measurements of grain emissivity improve, it will become necessary 
to assume specific models in order to estimate how different grain species 
are weighted in measurements of flux density and emissivity. The problem 
should be readily manageable for two-component models. 

As we will show in later sections (4, 5, 6), thermal emission mea- 
surements may themselves provide tests of alternative models. In particular, 
it may prove feasible to exploit the separation in temperature according 
to emissivity, the phenomenon at the crux of the complications for mixtures, 
to test the relationship of dielectric and magnetic properties of the 
principle species. 
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3 CLOUD MASS 

In determining a cloud mass from thermdl emission one need not know or 
assume values of either a or Q. One measures the visual extinction, A v , and 
the submillimetre optical depth, t( v ), in a calibration object with visible 
background stars. One assumes that the ratio of A y to hydrogen column den- 
sity is the same as that for the intercloud medium and that that ratio and 
the ratio A y /x(v) for the calibration object can be applied to other (denser) 
clouds. 

We shall write explicitly the relationship between dust mass and cloud 
mass; then the relationship of hydrogen column density to x(v); then, qual- 
itatively, the conditions favoring accurate results. 

3.1 Rela£co)i6lu.p oj cloud maM to duAt If one introduces a 

ratio, Mg/Mj, of gas mass to dust mass, [M g £ M( cloud)] then one may 
rewrite equation (4) to give the mass of a cloud: 

M c = [F(v)D 2 /B(v,T)] [(4/3) a/Q(v) ] p [M g /M d ] (6) 

Purcell (1969) has shown that if one assumes an extreme value 
e(0) = °°, for the static dielectric constant of grain material, one can 
use the integrated absorption to place a limit on the fraction of the 
total volume that is solid material. When this limit is combined with 
measurements of hydrogen column density, one may derive a lower limit 
(^60) for Mg/M^, 

In principle, one could use measurements of dust volume 
and hydrogen column density to derive a value for M^/M d . What we shall 
pursue here is the use of flux density and temperature measurements com- 
bined with extinction measurements to estimate the entire factor 

C = [(4/3) a/Q(v)] p [M /M d ] (7) 

appearing in equation (6). 

15 


3.2 Rzlatiomkip oj cloud m- 66 to optical depth. Using the ratio 


R = A v /E(B-V) of visual extinction, A y , to colour excess , E(B-V) (the 
excess over unreddened values of the difference between blue and visual 
magnitudes) and tie ratio N(H + H 2 ) /E ( B- V) of hydrogen column density 
to colourexcess as measured in the intercloud medium (Bohlin, Savage 
and Drake 1978) ,and assuming that the same numbers apply to a thin cloud, 
one obtains the column density per unit extinction, N(H + H 2 )/A y . 

The submillimetre optical depth is given by 

t(v) = F(v)|>e? BUT)]' 1 (x«l) (8) 

where F(v) is the flux density within a beam of angular radius e, . 

In a cloud thick enough to allow a measurement of t(v) but 
thin enough to allow a measurement of A y one may obtain a ratio A v / T (v) 
from which to derive column densities in any cloud for which t(v) is known 
[assuming A v / T (v) £ constant and N(H+H 2 )/A v £ constant]. Thus 

N(H+H 2 ) = T (v) [A v / t (v)] [N(H+H 2 )/E(B-V)] [1/R] (9) 

2 2 

The mass of the cloud within the area ire D covered by the beam is 

% 

M c = Ue^D 2 ] [N(H-H 2 )m HM ] 

where m^ is the mass of a hydrogen atom and y is the ratio of total gas 
mass to hydrogen mass 1.36). With equations 8 and 9, this becomes 



M c 

= [F(v)D 2 /B( v ,T)]C 

(10a) 

where 

C = 

[A v /t(v)] [N(H+H 2 )/E(B-v)] (1/R] ryi 

(10b) 

or 

C = 

[N(H+H 2 )/ T ( v )]m H M. 

(10c) 


This analysis has been presented by Keene, Hildebrand and Whitcomb (1982) 
using the data of Whitcomb et al. (1981) for the reflection nebula NGC 7023. 
The numerical values are given in Table 1. 
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3.3 Concloaf-on to dd^cuAAlcm oj cloud mas- 6. When applying the factor C, 


evaluated by means of equation 10, to the measurement of a dense cloud, 
one need not be concerned, as in the case of dust mass, with chosing a 
frequency at which each grain species emits in proportion to its volume. 

On the other hand, one does not entirely skirt the problem of differences 
in temperature among different grain species (section 2.5). It is 

still necessary that the various species emit in the same proportions as 
in the calibration object. To satisfy that condition, the temperature 
and ultraviolet exposure of the calibration object should be as near as 
possible to those of the clouds to be investigated. It would be desirable 
to extend the calibration to visually opaque clouds by measuring the ex- 
tinction at, say, 2 urn and applying the known relationship of 2 ym to 
visual extinction. 

In choosing a frequency at which to measure and apply the calibration 
one seeks to minimize the dependence of flux density on temperature by 
minimizing hv/kT. To good approximation maps made at A % 350 ym follow 
the distribution in column density. Hot regions in large molecular clouds 
generally have relatively low mass and have little influence on submilli- 
metre emission. For example, a map of the SgrA molecular cloud at 
540 ym shows no maximum at the position of the galactic center (Hildebrand 
et al . 1978) . 

The column density mapped by the submillimetre flux density is, of course, 
the duAt column density. Whether the ratio of gas to dust is in fact constant, 
as implied by the calibration procedure, is a matter to be explored by very 
careful comparisons from cloud to cloud and point to point of thermal emission 
and molecular line emission taking into account all the factors other than 
gas column density which may influence the intensity of the observed lines 
(Frerking, Langer, and Wilson 1982). 
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4 DIELECTRIC CONSTANTS 


It is evident, (equation 5) that a determination of a and Q imposes a 
constraint on the mean value of the imaginary part of the expression 
[e(v)-l ]/[e(v)+2] where e(v) is the complex dielectric constant. 

If the extinction curve is known over a sufficient range of wave- 
lengths, one may also determine the value of [e(0)-l ]/[e(0)+2] where 
e(0) is the static dielectric constant. Purcell (1969) has shown that 
the static electric susceptibility x'(0) of a cloud thinly populated by 
grains is given by 


x'(0) = ( 2/ 7t ) Q f [x"(v)/v] dv (11) 

where x"(v), the imaginary part of the susceptibility, is related to the 
total cross section by 

p 

4Trx" = nXctj. = n(c/2uv)Tra ^ext^ v ^' 


Here, n = number of grains per unit volume, = total cross section, 
and Q ext = extinction efficiency (= emissivity, Q, at long wavelengths). 
For spherical grains 


X* (0 ) = n a [e(0)-l]/[e(0)+2) 


( 12 ) 


Combining these expressions we have 


(4n^a/c){ [e(0)-l ]/[e(0)+2]} = /" «' 2 Q ext ( v) dv. 03) 

(See Aanestad and Purcell 1973.) 
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The modification of (12) for nonspherical grains (Purcell 1969) is 
significant only if e(0) is found to be much larger than the values 
(< 4) for typical insulating materials and if the grains are very far 
from spherical (ratio of axes >4 or £1/4). 

The evaluation of the integral is now limited primarily by the 
uncertainty in Q ext (v) at infrared and submillimetre wavelengths. 
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5 EMISSI VITY 


In the preceding discussion we have emphasized the importance of 
the emissivity, Q(v), in estimating dust masses of clouds and dielectric 
characteristics of dust grains. We now review the principles involved 
in measuring Q(v) at long wavelengths. 


For (UV) wavelengths comparable to grain radii, the absorption cross 

section, o abs , of a grain is nearly equal to the geometrical cross section, 
2 

Ogeom = vCL » anc l> by Babi net's principle, the scattering cross section, 
°scat’ must a1so be near ^ e( l ua 1 to °geom' hence f° r these wavelengths 
the total cross section. 


a tot °abs + a scat’ 


(14) 


becomes 


°tot (UV) 4 2o geon.. 


(A % a), 


(15) 


[A resonance near 0.2 ym increases a. . (UV) to %3 a ]• 

tot geom. 

For wavelengths much greater than grain radii, 

°scat * 0 amJ 


o t 0 t (lR) * ° ah 5 (lR), >> »)• 


abs 


(16) 


If one divides these expressions through by a , one obtains the 

geom 

corresponding expressions for the extinction, scattering, and absorption 
efficiencies: 


and 


^ext ^scat + ^abs ’ 


Q e xt^ UV ) ^ ^ (° r somewhat more) ( A £ a) 


Q ext (lR) = Q =»k c (ir) 


abs 


(a >> a ) 


(17) 

(18) 
(19) 
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where, by Kirchhoff's law, Q abs (IR) = emissivity, Q(v). 

The determination of Q(v) .it X » a consists in measuring the 
ratio Q(v)/Q ext (UV) , and multiplying by the known value of Q ex t(UV). 

The measurement of this ratio, as discussed in detail by Whitcomb 
et al . (1981), can be carried out by two methods. Briefly they are as 
follows: In the first, one finis Q(v)/Q ext (UV) from the equivalent ratio 

of optical depths t(IR)/t(UV). The long wavelength optical depth, t(IR), 
is the ratio (<< 1) of the observed brightness of a (thin) cloud to the 
brightness of a blackbocy of the same temperature (equation 8). The UV optical 
depth is obtained by measuring the visual or near IR extinction of background 
stars. In the second method, one carries out an energy balance for grains 
in a reflection nebula heated by a UV-emitting star of known spectral 
type. The energy absorted fromthe UV flux of the star is set equal to 
the observed thermal emission. If one assumes that the star-to-grai n dis- 
tance is equal to the projected distance and neglects UV absorption between 
the star and the grain one obtains an upper limit to the energy absorbed and 
hence a lower limit to Q(v). 

In carrying out the energy balance one must use a value for the frac- 
tion of the incident UV energy absorbed by the grain. We have implicitely 
taken that fraction, Q abs ( uv )/Q ex t( uv )' to be in ec l ua tion 18. Whitcomb 
has noticed that the UV and IR observations of a reflection nebula can be 
combined to give a measure of the albedo, y = Q sca t/Q ex t> anc * hence °f 
(1 - y) = Q abs /Q ext - The probability, P, that a UV photon from an imbedded 
star will escape from a cloud is just the ratio of the UV luminosity to the 
total luminosity. 
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P = L(UV)/[L(UV) + L( IR) ] • 


( 20 ) 


But P is also given by 

P = Y <m> (21) 

where <m> is the mean number of times a photon is scattered before it 
escapes. If the mean scattering angle is small, as indicated by the obser- 
vations of Lillie and Witt (1976), then <m> * x. Using the measured 
values of P (equation 20) and x (equation 8), one obtains a value of y 
(0.36) in good agreement with that obtained from observations of the dif- 
fuse galactic light (0.40; Lillie and Witt 1976). 

In measuring an "average" Q(v), values for individual grains are 
weighted according to size as discussed in section 2.2. The problem of 
determining Q(v) in mixtures of grains is discussed in section 2.5. 


22 


6 GRAIN COMPOSITION: POLARIZATION OF EMITTED RADIATION 


We have thus far considered only average properties of grains. We 
have discussed individual grain species only to the extent of arguing 
that the mix of species is not important with regard to determinations 
of dust volume. With regard to the evolution of the interstellar medium, 
however, the question is fundamental. 

Models of dust grain composition are based largely on examination of 
the extinction curve and on indirect evidence from studies of the chemical 
composition of interstellar gas as discussed, for example, by Field (1973). 

In the model of Mathis, Rumpl , and Nordsieck (1977), graphite grains, sug- 
gested by an absorption feature at 0.22 pm (Hoyle and Wickramasinghe , 1962), 
and silicate grains, suggested by a feature at 10 pm (Gillett, Low, and 
Stein, 1968; Woolf and Ney, 1969; Gillman, 1 969) , are present in approximately 
equal numbers and account for essentially the entire grain population. If 
that model is correct, then the graphite grains will be at a higher temper- 
ature than the silicate grains since graphite is a better ultraviolet and 
visible absorber and a poorer far-IR and submillimetre emitter (Mezger, 

Mathis, and Panagia, 1982). Consequently the flux on the short wavelength 
side of the composite emission peak should be mostly from graphite and the 
flux on the long wavelength side should be mostly from silicates. 

The separation in temperature provided by the difference in the emis- 
siviities of the grains provides an opportunity to test the model by examining 
another distinguishing characteristic, the polarization properties. The 
0.22 pm feature (graphite?) shows no enhanced polarization (Gehrels, 1974). 
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The vicinity of the 10 pm feature (silicates?) shows strong polarization 
in some objects, probably due to grain alignment (Dyck and Beichman 1974). 
One should therefore expect to find stronger polarization on the long 
wavelength side of the emission peak than on the short wavelength side 
if the graphite-silicate model is correct. The polarization should be 
normal to that at 10 pin if it is correct that the 10 ym polarization is 
due to selective absoprtion by aligned grains. 

For other models one would expect qualitatively different results. 

If, for example, the carbon is amorphous rather than crystalline as pro- 
posed by Draine (1981) and by Czyzak, Hirth and Tabak (1982), then the 
carbon emissivity at X £ 100 ym would be higher than that of silicates. 

An amorphous carbon-silicate mixture would therefore exhibit stronger 
polarization on the short wavelength side of the emission peak. 

The results of the Cornell group (Gull et aK , 1980) and the 

University College, London group (Cudlip et^ , 1982) in their infra- 
red polarimetry of Orion are consistent with the interpretation that 
the short wavelength side of the spectrum is dominated by emission 
from grains of low magnetic susceptibility (e.g., graphite). The 
polarization reported by Cudlip ert al^. (1982) at A = 77 pm is 
2.2 ± 0.4 percent, nearly orthogonal to the polarization in the near- 
infrared. The question to be resolved is whether the polarization at 
A » 100 pm is larger or smaller than this value. 

It has long been recognized that clouds containing aligned grains 
should emit polarized radiation (e.g., Stein, 1966). It is now possi- 
ble that submillimetre polarimetry could become a practical technique 
for investigation of grain characteristics. 
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7 NUMERICAL VALUES OF SUBMILLIMETRE EMISSION PARAMETERS 


In Table 1 we present numerical values for the quantities which 
enter into estimates of dust mass, cloud mass, and grain characteristics 
as derived from submillimetre emission. To avoid any implication of final- 
ity, our values are labeled “1983 Chicago assumptions." There is no basis 
for assigning statistical errors. Instead, we give a range of values 
consistent with various assumptions and measurements. We do not attempt 
to evaluate the inaccuracy in the implicit assumption that values for 
one cloud can be applied to another. 

Nowhere have we assumed a value for the ratio, M^/M^, of gas mass to 
dust mass. The value one would infer from the coefficients given in the 
last two lines of the table (10/0.1 = 100) is subject to the large errors 
claimed for those coefficients. 

An important goal of this paper is to encourage improved measurements 
of the quantities presented in the table. 
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APPENDIX 


AueAgge pAo/ec-tccn a. convex object onto a plan e . 

S. Chandrasekhar arid R. G. Sachs have called my attention to a 
theorem by K. Schwarzchild which states that the average shadow of a 
convex body is just one quarter of the total surface area. 

The theorem was re-discovered by M. Schwarzschi Id during World War II 
while he was working with Chandrasekhar and Sachs at the army ballistics 
laboratory in Aberdeen, Maryland. Subsequently, Schwarzschild learned 
from G. Birkhoff of his father's earlier work (M. Schwarzschild, private 
communication). Sachs recalls that A. Wintner was the first to recognize that 
the theorem was a re-discovery. I have been unable, even with the help 
of these distinguished colleagues, to find the original publication. It 
can doubtless be found somewhere in the Goettingen Archives where Martin 
Schwarzschild has sent his father's papers. 

The theorem is easy to prove when one knows that there is a simple 
answer. I offer the following: 

Divide the convex surface into n elements with n/2 exposed, on the 
average, to the observer. Let a i = outward vector representing the ith element, 
and r = unit vector toward the observer. 

Then t.-r = projection of the ith element on the sky. 

The average value of that projection for all directions of a. in the 
hemisphere facing the observer, i.e. for all rotations of the grain for 
which the ith element is visible, is just a./2*. Hence the mean projected 
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area, a, of the (average) n/2 visible elements is 


n/2 n 

a = £ a./2 = (1/4) £ a H 

i=l 1 i =1 1 


or 

o = A/4 


where A = total surface area. 


* a^ cose (2iTsine) dd/ Q f v ^ (2Trsine) de = a^/2. 
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Table I 

SubmilLime £> te Emtb^^ion F actosu 


Factor 


Sect, or 
eqn. in text 

1983 Chicago 
assumptions 


Range of values; 
comments 

6 

index giving fre- 
quency dependence of Q 
[F(v)«c/b(v,T)] 

§2 

1 for 50 ym <\<_ 

2 for X> 

250 ym 
250 ym 

^1 at 100 ym to 
_> 2 at 1000 ym 

P 

density of grain 
material 


3 gm/cm^ 


amorphous C: 1.9 
graphite 2.25 

enstatite 3.1 

Si C 3.2 

olivine 3. 2-3.4 

magnetite 5.0 

a 

ave. grain radius^ 
weighted by n(a)a 

§2.2 

0.1 ym 


0.25 ym-0.01 ym 
+smaller sizes 

Q(UV) 

grain emissivity at 
0.15 ym to 0.30 ym 

§5 

3 


2. 5-3.5 

Q(uv) 

Q 1 25 

ratio of Q(UV) to grain 
emissivity at 125 ym 

§5 

4000 


> 700 
1250 

700-2500 if 6 = 2* 
2000-5000 i f e = 1* 
* at * > 50 um 

Q 250 

grain emissivity at 250 
(assuming 6 = 1) 

ylTl 

( 125/250) 1 (3/4000) 
= 1/2666 



CONJTJMUE.U ONJ bJEKT PA G-EL 


N(H+H 2 ) 

T x 

ratio of hydrogen 
column dens, to opt. 
depth at wavelength \ 

eqn. 9 

1 . 2 x 10 ^ atoms/cm^ 
at 400 ym ? 

(1.2xl0 25 )(250/400r 

atoms/cm^ at 250 ym 

24 2 

7-26x10 atoms/cm 

at 400 ym; 

depends on 6 

C 250 

coeff. for estimates 
of cloud mass (250 pm) 
[C =C 25 Q (x/ 250) b ] 

eqn. 10 

2 

10 gm/cm 

value probably good 
within a factor 2 

( 4 / 3 ) ap 
Q 250 

coeff. for estimates 
of dust mass (250 ym) 
[Q =Q 250 (250/X)B] 

eqn. 4 

2 

0.1 gm/cm 

value probably good 
within a factor 3 or 
4 for measured cloud 
(NGC 7023) 


CJ 

o 


a. Schwartz 1982; Eri ckson ejt aj_. 1981 

b. Mathis, Rumpl , and Nordsieck 1977; Mathis and Wallenhorst 1981 

c. Lillie and Witt 1976 

d. Aanestad 1975 

e. Keene 1981 

f. Harvey, Thronson, and Gatley (1980): result converted to correspond to 125 pm figure of Whitcomb et al_. (1981) 

g. Whitcomb, Gatley, Hildebrand, Keene, Sellgren, and Werner 1981 

h. Keene, Hildebrand and Whitcomb 1982 

i. Derived from expression for N(H+H2)250’ c ^* e Q ua ti° ns 10 and 11 

j. From assumed values of a, p and Q 2 gQ. 


REFERENCES 


Aanestad , P. A., 1975. Absorptive properties of silicate core-mantel 
grains, Astrophys. J. 200 , 30. 

Aanestad, P. A. & Purcell, E. M. , 1973. Interstellar grains, 

Ann. Rev. Astron. Astrophys. 1_1_> 309. 

Aanestad, P. A. & Kenyon, S. J., 1979. Temperature fluctuations and the 
size distribution of interstellar dust grains, Astrophys. J. 230 , 

771. 

y 

Andriesse, D. C., 1977. Radiating cosmic dust. Vistas in Astronomy 2J_, 107. 

Bohlin, R. C. , Savage, B. D. & Drake, J. F. , 1978. A survey of inter- 
stellar HI from L absorption measurements. II. Astrophys. J. 224 , 132. 

Carrasco, L., Strom, S. E. & Strom, K. M., 1973. Interstellar dust in 
the Rho Ophiuchi dark cloud, Astrophys. J. 182 , 95. 

Coyne, G. V., Gehrels, T. & Serkowski, K. , 1974. Wavelength dependence of 
polarization. XXVI. The wavelength of maximum polarization as a 
characteristic parameter of interstellar grains, Astron. J. 79_ , 581. 

Cudlip, W., Fruniss, I., Kind, K. J., & Jennings, R. E., 1982. Far-infrared 
polarimetry of W51A and M42, Mon. Not. R. astr. Soo. 200 , 1169. 

Cyzak, S. J., Hirth, J. P. & Tabak, R. C. , 1982. The formation and proper- 
ties of grains in the interstellar medium. Vistas in Astronomy 25^, 337. 

Drain, B. T. , 1981. Infrared emission from dust in shocked gas, 

AAtnophy* . J. 245 , 880 . 

Dyck, H. M. & Beichman, C. A., 1974. Observations of infrared polari- 
zation in the Orion Nebula, Aifriophy*. J. 194, 57. 

Erickson, E. F., Knacke, R. F., Tokunaga, A. T. & Haas, M. R., 1981, 

The 45 micron H 2 0 ice band in the Kleinman-Low nebula, 

AAtAopky-b . J. 245 , 148. 

Field, G. B., 1973. The composition of interstellar dust, in The Dusty Universe , 
eds. Field, G. B. & Cameron, A. G. W., Neal Watson Academic Publications, 

New York, 89-112. 


31 


Frerking, M. A., Langer, W. D. , & Wilson, R. W., 1982. The relationship 
between carbon monoxide abundance and visual extinction in inter- 
stellar clouds. Astrophys. J. 262 , 590. 

Gehrels, T., 1974. Wavelength dependence of polarization. XXVII. Inter- 
stellar polarization from 0.22 to 2.2/mm. Astrophys. J. 79, 590. 

Gillett, F. C., Low, F. J., & Stein, W. A., 1968. Stellar spectra from 
2.8 to 14 microns. Astrophys. J. 1 54 , 677. 

Gilman, R. C., 1969. On the composition of circumstellar grains. 

Astrophys. J. ( Letters ) 1 55 , LI 85. 

Greenberg, J. M. & Hong, S. S., 1973. Effects of particle shape on 
volume and mass estimates of interstellar grains, in The Dusty 
Universe , eds., Field, G. B. & Cameron, A. G. W., Neal Watson 
Academic Publications, Inc., New York, 131-154. 

Gull, G. E., Russell, R. W., Melnick, G., & Harwit, M. , 1980. Far-infrared 
polarization of the Kleinmann-Low nebula. Astron. J. 85, 1979. 

Harvey, P. M., Thronson, H. A., Jr., & Gatley, I., 1980. A far-infrared 
study of the reflection nebula NGC 2023. Astrophys. J. 235 , 894. 

Hildebrand, R. H., Whitcomb, S. E., Winston, R., Stiening, R. F., 

Harper, D. A., & Moseley, S. H., 1977. Submillimeter photometry 
of extragalacti c objects. Astrophys. J. 21 6 , 698. 

Hildebrand, R. H., Whitcomb, S. E., Winston, R. Stiening, R. F., 

Harper, D. A. , & Moseley, S. H., 1978. Submillimeter observations 
of the galactic center. Astrophys. J. (Letters) 219 , L 1 0 1 . 

Hong, S. S. & Greenberg, J. M. , 1980. A unified model of interstellar 
grains: a connection between alignment efficiency, grain model 

size, and cosmic abundance. Astron. Astrophys. 88, 194. 

Hoyle, F. & Wickramasinghe, N. C. , 1962. On graphite particles as inter- 
stellar grains. Mon. Not. R. Astr. Soo. 124 , 417. 


32 


Jackson, J. D. , 1961. Introduction to dispersion relation techniques 

In Dispersion Relations, ed. , Screaton , 6. R. Interscience Publishers, 
Inc., New York, Chapter 1. 

Keene, J., 1981. Fir infrared observations of globules. Astrophys. J. 

245 , 115. 

Keene, J., Hildebrand, R. H. & Whitcomb, S. E. , 1982. A high resolution 
map of OMC-1, AttAophy*. J. (Le£tm), 252, LI 1 . 

Leung, C. M. , 1975. Radiation transport in dense interstellar dust clouds. 

I. Grain temperature, AitAophy* . J. 199 , 340. 

Lillie, C. F. & Witt, A. N., 1976. Ultraviolet photometry from the 

orbiting astrophysi cal observatory. XXV. Diffuse galactic light 
in the 1400-4200 A region and the scattering properties of inter- 
stellar dust grains, AbtAopky*. J. 208 , 64. 

Mathis, J. S., Rumpl , W. & Nordsieck, K. H. , 1977. The size distribution 
of interstellar grains, AbtAophy*. J. 217 , 425. 

Mathis, J. S. & Wallenhorst, S. G. , 1981. The size distribution of 
interstellar particles. III. Peculiar extinctions and normal 
infrared extinction, A itAophy*. J. 244, 483. 

Merrill, K. M. , Russell, R. W. & Soifer, B. T., 1976. Infrared observations 
of ices and silicates in molecular clouds, AitAophy*. J. Z07_, 763. 

Mezger, P. G. , Mathis, J. S. & Panagia, N., 1982. The origin of the 

diffuse galactic far infrared and submillimeter emission, A itAon. and 
At>tAopky & . 105 , 372. 

Purcell, E. M. , 1969. On the absorption and emission of light by 
interstellar grains, A5 tAophyi>. J. 158 , 433. 

Purcell, E. M. , 1976. Temperature fluctuations in very small interstellar 
grains, A& tAophy*. J. 206 , 685. 

Schwartz, P. R. , 1982. The spectral dependence of dust emissivity at 

millimeter wavelengths, A ttAophy*. J. 252 , 589. 

33 


Serkowski, K. , Mathewson, D. S. & Ford, V. L., 1975. Wavelength dependence 
of interstellar polarization and ratio of total to selective extinction, 
AAtn.oph.y6. J. 196, 261. 

Stein, W. , 1966. Infrared radiation from interstellar grains, 

A 6tsiophyA. J. 1 44 , 318. 

Whitcomb, S. E., Gatley, I., Hildebrand, R. H., Keene, J., Sellgren, K. 

& Werner, M. W., 1981. Far infrared properties of dust in the reflec- 
tion nebula NGC 7023, AAfrwpky*. J. 246 , 416. 

Witt, A. N., 1979. Ultraviolet scattering properties and the size 

distribution of interstellar grains, fatnopky*. and Space See. 65^ 21. 

Woolf, N. J. & Ney, E. P., 1969. Circumstellar infrared emission from 
cool stars. Astrophys. J. (Letters) 155 , L181 


34 


1. Report No. 2. Government Accession No. 

NASA TM-84356 

3. Recipient's Catalog No. 

4. Title and Subtitle 

THE DETERMINATION OF CLOUD MASSES AND DUST CHAR- 
ACTERISTICS FROM SUBMILLIMETRE THERMAL EMISSION 

5. Report Oate 

6. Performing Organization Code 

7. Author(s) 

Roger H. Hildebrand 

8. Performing Organization Report No. 

A-9313 

10. Work Unit No. 

352-02-03 

9. Performing Organization Name and Address Department of Astronomy and 
Astrophysics, Department of Physics, The Enrico 
Fermi Institute, The University of Chicago, 

Chicago, Illinois 60637 

11. Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Memorandum 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

14. Sponsoring Agency Code 


15. Supplementary Notes Preprint Series # 9. Supported by NASA grants. Point of 

Contact: L. C. Haughney, Ames Research Center, M/S 211-12, Moffett Field, 

Calif. 94035, (415) 965-5339, FTS 448-5339. 


16 Absract 


The purpose of this paper is to review the principles by which the 
dust masses and total masses of interstellar clouds and certain character- 
istics of interstellar dust grains can be derived from observations of 
far-infrared and submillimetre thermal emission. To the extent possible, 
the discussion will be independent of particular grain models. 


17. Key Words (Suggested by Author(sf) 

Astronomy 

Interstellar matter 
Interstellar abundances 
Far-infrared and submillimetre 

18. Distribution Statement 

Unlimited 

STAR Category - 89 


19. Security Oassif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* 

Unclassified 

Unclassified 

36 

A03 


'For sale by the National Technical Information Service, Springfield, Virginia 22161 


